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Abstract 

We derive the LO DGLAP evolution equation for the full Mellin moments of the truncated at xo first 
moment of the nonsinglet parton distribution. This "moment of moment" approach allows to determine the 
small-xo behaviour of the truncated first moment. We compare our predictions to results obtained from 
a;— space solutions for parton distributions with use of the Chebyshev polynomial method and to solutions of 
the evolution equations for the truncated moments proposed by other authors. The comparison is performed 
for different input parametrisations for 10 -5 < xo < 0.1 and 1 < Q 2 < 100 GeV 2 . We give an example of 
application to the determination of the contribution to the Bjorken Sum Rule. 
PACS 12.38.Bx Perturbative calculations, 11.55.Hx Sum rules 

1 Introduction 

Deep inelastic scattering (D1S) experiments provides knowledge about the internal structure of the nucleon. 
Measurements on proton, deuteron and neutron targets allow also verification of sum rules e.g. ^ 12 - 
basic relations in QCD. Sum rules for unpolarised and polarised structure functions refer to their Mellin 
moments. Particularly important role in the experimental and theoretical QCD tests play first moments 
of the parton densities, which have physical interpretation and can be determined from the data. From an 
experimental point of view, however, the accurate verification of sum rules is unreliable. Determination of 
the sum rules requires knowledge of the structure functions over the entire region of the Bjorken variable 
x 6 (0; 1). The lowest limit of x in present experiments is about 10 and the limit x — > 0, which means 
that the invariant energy W 2 of the inelastic lepton-hadron scattering becomes infinite (W 2 = Q 2 (l/x — 1)), 
will never be attained. Therefore it is very hopeful in the theoretical analysis to deal with truncated instead 
of with full moments of the structure functions. This enables to avoid uncertainties from the unmeasurable 
x — > region. The most familiar theoretical approach, which describes scaling violations of parton densities 
in perturbative QCD has been formulated by Dokshitzer, Gribov, Lipatov, Altarelli and Parisi (DGLAP) 
0. Evolut ion equations for the truncated at low xo moments of parton distributions are however more 
complicated than in a case of the full moments. These are not diagonal and each n-th truncated moment 
couples to (n + j)-th (j > 0) truncated moments g|. For n > 2 the series of couplings to higher moments is 
very fast convergent. Even for small (m = 4) number of terms in the expansion of the truncated counterpart 
of the anomalous dimension G„ , the higher moments can be calculated with excellent accuracy. First moment 
is more sensitive to the truncated point xo and the convergence of G n for n — 1 is weaker than for the higher 
moments. Nevertheless, it has been shown in |5], that for more terms of the G„ expansion (m ~ 30), the 
uncertainty in the determination of the first moment at Xo < 0.1 and Q 2 = 10 GeV 2 does not exceed 5% 
independently on the input parametrisation. However, this increase with m of the accuracy does not proceed 
infinitely. Numerical errors, which occur for larger m (dependently on xq) make further improvement of the 
precision impossible. It may be very useful to discuss methods of theoretical determination of truncated first 
moments of parton distributions because these predictions would be directly verified experimentally. In this 
paper we present an approach, in which we compute the truncated at small xo first moments using the inverse 
Mellin transform of their full moments. In other words, our method is based on the solutions for the full n-th 
moments of the truncated first moment of the parton distribution. This "moment of moment" technique 
would be a complementary one to other known methods in determination of sum rule contributions. 

The content of this paper is as follows. In Section |5]we recall the ways of computing the truncated 
moments within DGLAP approximation. Thus Sect. |2~T1 contains a brief description of the Chebyshev- 
polynomial approach for x-sp&ce solutions of DGLAP evolution equations. In Sect. 12.21 we recall the idea of 
DGLAP evolution equation for the truncated moments, which underlies our "moment of moment" technique. 
This modified method is presented in Sect. 12.31 For simplicity we consider the quark nonsinglet LO evolution. 
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In Section 0we compare solutions for the truncated at small xo first moment of the nonsinglet structure 
function, obtained with use of these mentioned earlier approaches. As an example we explore the Bjorken 
Sum Rule (BSR). In Section [I] we present predictions for the low- a: contribution to the BSR together with 
experimental constraints. Finally, in Section U^we summarise our results. 



2 Determination of the truncated Mellin moments of parton 
distributions. 

For (full) Mellin moments of parton distributions p(x,Q 2 ) 

1 

M[p]{n,Q 2 ) = / dxx n ~ x p{x,Q 2 ) (1) 



the DGLAP evolution equation can be solved analytically. This is because in the moment space n one obtains 
simple diagonalised differential equations. The only problem is the knowledge of the input parametrisation for 
the whole region < x < 1, what is necessary in the determination of the initial moments M [p] (n, Q 2 = Qo): 



i 



M[p](n,Q 2 )= / dxx^pfaQl). (2) 



o 



Using the truncated moments approach one can avoid the uncertainties from the region x — ► 0, which will 
never be attained experimentally. 

The truncated at xq Mellin moment of the parton distribution p(x,Q 2 ) is defined as 

i 

M[p](x ,n,Q 2 ) = I dxx^pix^Q 2 ). (3) 



From the theoretical point of view, there are two ways to avoid the problem of dealing with the unphysical 
region x — * 0. The first one is to work in a;-space and obtain directly the evolution of parton distributions (not 
of their moments). The most known methods for solving the Q 2 evolution equations for parton distributions 
in a>space are brute- force [Bj, Laguerre-polynomial |7j or Chebyshev-polynomial [S| approaches. In this way, 
the truncated moment can be simply found by integrating the x-space solutions p(x,Q 2 ) over the cut range 
xo < x < 1 (see Sect. 12. it . An alternative way is use the evolution equations directly for truncated moments. 
The appropriate DGLAP evolution equations for the truncated moments have been derived in Authors 
have shown that these equations, though not diagonal, can be solved with a quite good precision for n > 2, 
even for a very small number of terms in the expansion series. In a case of the first moment, the accuracy 
is worse and more terms in the G n expansion must be taken into account. We briefly recall the idea of 
solving the evolution equations for truncated moments in Sect. 12.21 Basing on this idea, we have derived 
the evolution equation for truncated first moment in a diagonal form. The appropriate integro-differential 
equation contains only one function - gi(x, Q 2 ), which denotes truncated at x first moment, without coupling 
to the other, higher moments. Then, using the full Mellin moments approach, we have found the small- 
x = xo behaviour of the function qi(x,Q 2 ), which is simply the truncated at low-xo first moment. Detailed 
description is given in Sect. 12.31 



2.1 LO DGLAP evolution equations with use of the Chebyshev- polyno- 
mial expansion. 

The Chebyshev polynomials technique |0] was successfully used by J.Kwiecihski in many QCD treatments 
e.g. [Sj. Using this method one obtains the system of linear differential equations instead of the original 
integro-differential ones. For example, in order to solve the LO DGLAP evolution equation for the nonsinglet 
parton distribution p = q NS : 



dp(x,t) _ a a {t) 
dt ~~ 2tt 



(J) 



one has to expand functions p(x, t) into the series of the Chebyshev polynomials. In this way, the integration 
over z in the evolution equation @ can be performed, what leads to the system of linear differential equations: 

*^ = 5>«p(*,,t). (5) 

3 = 



This system can be solved by using the standard Runge-Kutta method with initial conditions given by 
the input parametrisation p(xj,ta). Truncated at xo moments are simply computed numerically via the 
integrating The Chebyshev expansion provides a robust method of discretising a continuous problem. 
This allows computing the parton distributions for "not too singular" input parametrisation in the whole 
x 6 (0; 1) region. More detailed description of the Chebyshev polynomials method in the solving the QCD 
evolution equations is given e.g. in Appendix of 0. 

2.2 Evolution equations for truncated moments. 

LO DGLAP evolution equation for the truncated at xq Mellin moment of the nonsinglet structure function 
found by Forte, Magnea, Piccione and Ridolfi (FMPR) 2] nas a form: 

dM[p](x ,n,t) as (t) f n f Xo \ 

1 dyy p(y,t)Gn — , (6) 



at An J \ y 

where 



* = (7) 



and 

i 



^QCD 



£?»(*) = J dyy"- L P qq (y) (8) 

Z 

is the truncated anomalous dimension. Expanding G n {xo/y) into Taylor series around y = 1 gives 

= ^ £ ^ <*0 M [p] (xo, » + i, (9) 

3=0 

Eq.@ is not diagonal but each n-th moment couples only with (n + j)-ih (j > 0) moments. Because 
the series of couplings to higher moments is convergent and furthermore the value of (n + j)-th moments 
decreases rapidly in comparison to the n-th moment, one can retain from @ the closed system of m + 1 
equations: 

dM[p](x ,n,t) _ as (t) 
dt ~ 2tt 

N +m 

x E D^ +m - n \x )M[p](x ,j,t). (10) 

j — n 

Here 

N <n<N +rn, (11) 
where iVo is the lowest considered moment and D is a triangular matrix. The solution of 1101 has the form: 

N +m 

M[p](x ,n, t) = A nk (xo)M[p](x ,k,t) 

fc = n + l 
N +m 

+ (M[p](x ,n,to)- J2 A nk (x )M[p](x ,k,t )) 

k = n + l 

x eatp^c/D^^ioJln^) . (12) 

Matrix elements D^j(xo) and Aij(xo) are given in In results l|12fl have been compared to those, 

obtained with use of the Chebyshev-polynomial technique. The agreement of both approaches is excellent 
for higher moments (n > 2) and not too large xq < 0.1, even for a small number of terms (m = 4) in the 
truncated series 1121 . It also does not depend strongly on the scale Q 2 or the input parametrisation. In a 
case of the truncated first moment, a similar accuracy requires more terms (m > 30) taken into account. 
In the next section we estimate the small- xo behaviour of the truncated first moment M [pi (xo, n= l,t). 

2.3 Small-xo behaviour of the truncated first moment. 

We would like to present a possible way of an approximate determination the small-so behaviour of the 
truncated first moment Mumxo, Our starting point is the evolution equation ©, which for the first 
moment has the form 

l 

d qi (x ,t) _a s (t) r dyp{Vit)Gl (Eo\ (13) 



dt 2tt / ; V y 



Here we denote in short the truncated first moment by qi(xo,t): 

q 3 (x ,t) = M[p](x J,Q 2 ). (14) 
Inserting Gi (z) in the LO approximation: 

Gi(*) = fln(l-*) + | (z+i), (15) 



3 v J 3 \ 2 / 
we obtain 

dqx(x ,t) 2a 3 (t) 

7T, = —5 [xo qo(x Q ,t) 

.2 M „fc 



+ %_ 1 (x 0j t)-2^^ gi _ fc (xo,t)]. (16) 



it 



In the above formula we have used the expansion of 
ln(l-z): 

OO k 

hi-*) = -J2t- (17) 

fc=l 

Taking into account that 

1 

Qj(xo,t) = x J ~ 1 q 1 (x ,t) + (j - 1) J dyy^ 2 qi (y,t), (18) 

we find 

9gi(a;o,i) _ 2a s (t) 



OO 

^(-o,t)(|-2Et 



at 3tt l " v ' ^2 ^ k> 

fc=i 

1 

The above result is LO DGLAP evolution equation for the truncated first moment of the nonsinglet parton 
distribution. This formula contains only qi and there is no mixing between qi and higher or lower moments 
qj. Therefore we can adopt to 1191 well known approach and write the evolution equation of qi in the 
moment space, what reads as follows: 

dM\qi] (n,t) 2a (t) r n 

— = ^r H ^ M l^ ( n > *)• ( 2 °) 

H(n) is given by 

H ( n )=2Y(-L r -l) + l--± — . (21) 

y J Z-^\n + k kj 2 n + 1 n + 2 v ' 

k=l 

In this "moment of moment" approach (MM) we have obtained a simple equation for a some complicated 
structure MTgi] (n, i), which is the (full) n-th moment of the truncated first moment of the parton density, 
namely 

1 1 

M[qi]{n,t) = J dxx"- 1 J dyp(y,t). (22) 

X 

The solution of J20j is 



M [<?i] (n, t) = M [qi] (n, to) exp[a(t)H(n)], (23) 
where M \qi \ (n, to) is a initial value of M \qi\ for the low scale to: 

1 1 

M[qi](n,to) = J dxx"- 1 J dyp{y,t ) (24) 

X 

and 

^=33^^ 

with Nf - the number of the quark flavours. Hence qi(x,t) is given by the inverse Mellin transform 

c+ioo 

qi(x,t) = — / dnx- n M[ qi ](n,to)exp[a(t)H(n)]. (26) 



The behaviour of q\{x,t) for x — > is governed by the leading (i.e. rightmost) singularities of (n, to) 

as well as H(n) in the n complex plane |10|. If we parametrise the input parton distribution p(x,to) in the 
general form 

p(x,t )~x a i(l-x) a \ (27) 

we obtain 



M \qx] (n, t ) ~ - /3(ai + 1, o a + 1) - V" 



fc! 

fc=0 

r(o 2 + i) 



(28) 



T(a 2 + 1 - fc)(ai + 1 + k)(n + 01 + 1 + fc) ' 
r(z), j3{zi,Z2) in 1281 are Euler functions and fc max depends on a 2 in the following way: 

, f oo for fractional a 2 > 

maa: I «2 for whole a 2 > 

One can see from 12 U and 1281 that H(n) has the rightmost pole at n = — 1, while M [51] (n, to) at n = 
and n = — ai — 1. In this way, for the nonsingular at small- a; shape of the starting distribution p(x,to) 
(ai = 0), the simple pole at n = and the essential singularity at n = — 1 are the leading ones. Then the 
small-xo behaviour of the truncated first moment can be determined by the method of steepest descent. We 
find 

qi (x,t) » — — - .fEp(a2 + l,z) x z 1 - 5 [z + 2a(t)]-°- s 

a-2 + 1 V Z7T 

x exp (a(t)H(z - 1) + 0.5^1 - 4a(t) ln(x )j , (30) 

where 

l + y/l-4a(t)ln(x ) 

2ln(x ) ■ {6 > 

If we consider more singular input parametrisation p(x,to) (01 < 0), this singular small- x behaviour remains 
stable against LO Q 2 QCD evolution. In this case the approximate behaviour of the truncated first moment 
qi(x, t) is governed by the leading simple poles of M [51] (n, to), situated at n = and n = — ai — 1: 

qi(x,t) « /3(ai + l,a 2 + 1) X ' ' 



^ fc! 

fc=0 

[ ' ' " -ea;p[a(i)ff(-ai - 1 - fc)]. (32) 



T(a 2 + 1 - fc) (01 + 1 + k) 

In the next section we compare results 13011 - 132(1 with those, obtained within approaches, described in 12. II 
and 



3 Results for the truncated first moment within three differ- 
ent approaches. 

In this section we present numerical results for the truncated first moment of the nonsinglet parton distribu- 
tion. Predictions obtained with use of different methods are denoted via CHEB, FMPR or MM, according 
to the applied approach (see Sect . 12 . ll 12 . 21 respectively') . 
Thus qi HEB (xo,t) results from integrating 

1 

CHEB 1 ,\ I j CHEB, ,\ / 00 \ 

qi (xo,t) = / dxx p (x,t), (33) 

XQ 

where p CHEB (x, t) is the solution of Q-JSJl, while q FAIPR (xo,t) is implied by 1121 and has the form 

m + l 

q[ MPR (x ,t) = £ A lk (x )q k (x ,t) 
k=2 

/ m+l 

+ I qi(x ,t ) - A lk (x )qk(xo,to) 



x exp 



(c,# (*„)]*£), (34) 



where q k (xo,t) is defined by 114H . 

Analytical approximate solutions 130H - I32L describing the low-xo behaviour of qi(xo,t) within "moment 
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Figure 1: Small xo behaviour of the truncated at xo first moment of the nonsinglet spin structure function 
g^ s in a case of the flat input (|35Jl • A comparison of (|3C)|) (MM) with the predictions based on the Chebyshev 
polynomials method (CHEB) and FMPR approach l|34|) for two values of the number of terms in the truncated 
series: m — 4, m = 30 is shown. The Bjorken Sum Rule is normalized to 1. 



of moment" (qi IM ) approach are compared with q\ HBB and q FMPR results in Figs. 1-6. We consider two 
values of the number of terms in the truncated series for the FMPR-m predictions (m = 4, m = 30). We 
show the results for the truncated at xo first moment of <?i , which is a contribution to the Bjorken Sum 
Rule |2] ■ In presented figures the total Bjorken Sum Rule is normalized to 1. We use two different inputs at 
Qo — 1 GeV 2 , namely: 

p{x,Ql) =g^ s (x,Q 2 ) ~ (l-xf, (35) 
p(x,Q 2 )=g? s (x,Q 2 ) ~ x-° A {l-xf-\ (36) 

More singular at small- a; parametrisation 1361 incorporates the latest knowledge about the low- a; behaviour 
of the polarised structure functions Truncated at xq = 0.01 integral J dxg± s is reduced by about 

8% for the Regge input 13511 and about 20% for 13611 in comparison to the total Bjorken Sum Rule. Figures 
1-6 show that for small- xo < 0.01 there is a good agreement between the MM results 1301 - 1321 and the 
predictions, obtained with use of the Chebyshev polynomial approach, which can be regarded as reliable. 
The percentage error 



qUxo ,Q 2 )-g? HEB (x ,Q 2 ) 
q? HEB (xo,Q 2 ^ 



E a (x , Q ) = 7 " U ^" ■ 100%, (37) 



where a denotes MM, FMPR-4 or FMPR-30 results, is about 3% in the case of MM solutions for x = 0.01 
and Q 2 — 10 GeV 2 , independently on the shape of the input parametrisation. Similar accuracy give taking 
into account more terms (m = 30) in the truncated series within FMPR approach (FMPR-30), while for 
m = 4 the error 1371 is about 4% in the case of the flat input and 6% in the case of the more singular one. 
E a (x ,Q 2 ) grows slowly with increasing Q 2 (see Figs. 5-6) and for Q 2 = 100 GeV 2 we get: E MM (x = 
0.01, Q 2 ) w E FMPR ~ 30 (x = 0.01, Q 2 ) fa 4% for the flat parametrisation and 7% respectively for the more 
singular input. Note, that for the truncation points xo < 0.01 our approximate analytical solutions I30t - I32j 
are as reliable as the FMPR — 30 predictions and more exact than the FMPR — 4 results. This does not 
depend either on the shape of the input parametrisation nor on the value of Q 2 . 
In the next section we determine the low-x contribution to the Bjorken Sum Rule. 



4 Low- a: contribution to the Bjorken Sum Rule. 

Among all moments of structure functions, the Bjorken Sum Rule (BSR) |5] is one of the convenient tests 
of QCD. BSR is a fundamental relation for polarised scattering, describing a relationship between spin 
dependent DIS and the weak coupling constant defined in neutron /3-decay. In the limit of the infinite 
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Figure 2: Small xq behaviour of the truncated at xo first moment of the nonsinglet spin structure function g^ s 
in a case of the steep input l|36|) . A comparison of ilML'l) (MM) with the predictions based on the Chebyshev 
polynomials method (CHEB) and FMPR approach 1)34(1 for two values of the number of terms in the truncated 
series: m = 4, m = 30 is shown. The Bjorken Sum Rule is normalized to 1. 
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Figure 3: Truncated first moment qf IM (|32fl as a function of ai in the input parametrisation p(x, Qq) ~ 
x ai (l — x) 3 for fixed xq = 0.01 and xq = 0.05. A comparison with the predictions based on the Chebyshev 
polynomials method is shown. 




Figure 4: Truncated first moment qf IM (|3U[) as a function of a-i in the input parametrisationp(a;, Qo) ~ (1 — x ) a2 
for fixed xq — O.Of and xq = 0.1. A comparison with the predictions based on the Chebyshev polynomials 
method is shown. 
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Figure 5: Q 2 dependence of the truncated first moment q± M (|30[l at fixed xo = 0.01 for input parametrisation 
l|35|) . A comparison to the predictions based on the Chebyshev polynomials method and FMPR approach l|34|) 
for two values of the number of terms in the truncated series: m = 4, m = 30 is shown. 
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Figure 6: Q 2 dependence of the truncated first moment q^ M l)32|) at fixed xq — 0.01 for input parametrisation 
1(36(1 . A comparison to the predictions based on the Chebyshev polynomials method and FMPR approach 1)34(1 
for two values of the number of terms in the truncated scries: m = 4, m = 30 is shown. 



momentum transfer Q , the BSR reads: 



r?-r? = J dx (#(*)- 



6 g v ' 



(38) 



where gv and gA are the vector and axial vector couplings. From recent measurements |12| gA/gv = 
1.2695 ± 0.0029. BSR refers to the first moment of the nonsinglet spin dependent structure function g\ S : 

g? s (x, Q 2 ) = 9 {{x, Q 2 ) - £(x, Q 2 ), (39) 

where g\ and g" are spin structure functions for proton and neutron. The asymptotic relation 13811 at finite 
Q 2 S> J^qcd takes a form with pQCD corrections: 



dxgi S (x,Q 2 ) = \ — 
b gv 



1 



— - 3.583 

7T 



(40) 



The validity of the sum rule is confirmed in polarised DIS at the level of 10% |13l 1141 . Ewaluation of the 
sum rules requires knowledge of polarised structure functions over the entire region of x: < x < 1. The 
experimentally accessible x range for the spin dependent DIS is however limited (0.7 > x > 0.003 for SMC 
data |13| . 0.6 > x > 0.023 for HERMES data and therefore one should extrapolate results to x — 
and x = 1. The extrapolation to x — > 0, where structure functions grow strongly, is much more important 
than that to x — > 1, where structure functions vanish. The extrapolation towards x — suffers from large 
uncertainties, being essentially dependent on the used QCD fit. "Flexibility" of the chosen parametrisation 
appears in the agreement with the experimental data, giving however enough freedom in the unmeasured 
regions |16|. In a case of the BSR this allows for a significant reduction of the low- a; contribution. 
In our approach we can test how the small- a; contribution to the BSR depends on the different (less or more 
steep) input parametrisations at the initial scale Q$ = GeV 2 : 



NS 1 

9i (x, 



) =Tj X ai (1 — x)" 



(41) 



Here r] is a normalization factor: 
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Figure 7: Low- a; contribution to the Bjorken Sum Rule 143(1 for different a\ and 02 in the input parametrisation 
a; = 0.01 and Q 2 = 5GeV 2 . 



The exponent ai controls the behaviour of the structure function g^ s as x — * and the factor (1 — a;)™ 2 
ensures the vanishing of gi S at x —> 1. The percentage contribution to the BSR, coming from small- x region 
< x < xo is defined as: 

Jdxg? s (x,Q 2 ) 

r(x , Q 2 ) = \ 100% (43) 

/ dx S f s (x,g2) 


The ratio r for xo = 0.01 varies from a few to tens percents for different configurations of — 1 < ai < and 
< a,2 < 8. The r-distribution at xo = 0.01 is shown in Fig. 7. One can see that the small- x contribution 
to the BSR grows with increasing 0,2 and decreasing a±. For the Regge flat parametrisation 1351 r is at the 
level of 5— 10%, what is significantly different from the result based on the input 1361 . where r ~ 20%. From 
theoretical analyses it is known that the small- a; behaviour of the nonsinglet polarised structure function 
gi S is governed by the double logarithmic terms i.e. (a s ln 2 x) n [T7|[TH1[TT]- This leads to the singular at 
low- a; form of gi S : 

9i S (x,Q 2 ) ~ x~ x (44) 

with A ~ 0.4. LO DGLAP approach with use of the singular input 1361 pretends the double logarithmic ln 2 x 
resummation. Low- a; experimental data |13l I15| clearly confirm the rise of g^ s in this region. However, the 
errors on the present data are too large to reliably support or contradict this x~ 0A behaviour. "Freedom" 
in the initial parametrisation to satisfy the small- a; experimental extrapolation of the BSR is seen in Fig. 8. 
The experimental data can be satisfactorily reproduced by e.g. nonsingular as x — > input ~ (1 — a:) 6 and 
by e.g. the singular one ~ a; _0,5 (l — x) 1 as well. SMC and SLAC measurements imply that 10 — 20% of the 
BSR comes from x values less than 0.01 I13II14|[T^| . Also HERMES data enable to determine the low- a; 
contribution to the BSR at 0.023 between 20 — 40%. Wide range of these estimations could be restricted by 
new spin data concerning this mystery and interesting small- a; region. 

5 Conclusions. 

In this paper we have compared results for the truncated at xo first moment qi of the parton distribution 
obtained within different approaches. Thus we have solved numerically LO DGLAP evolution equation for 
the nonsinglet function in the a>space using Chebyshev polynomial expansion and then after integrating over 
x we have got the prediction qi HEB , which can be treated as an exact one. Next, using evolution equations 
written in the moment space, we have found the closed system of m + 1 solutions for truncated moments, 
where obtained q% ~ m result is expressed by values of the m higher moments. Considered number of 
terms in the truncated series was m — 4 and m = 30. Working in the moment space we have also found 
an alternative way to determine the small- xo behaviour of the truncated first moment. Taking into account 
the relation between n-th and j-th truncated moment, we were able to derive the evolution equation, which 
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Figure 8: Constraints on the parametrisation of gf , implied by the experimental estimations of the low- a; 
contribution to the Bjorken Sum Rule. SMC, SLAC: 10% < r{x = 0.01,5) < 20% (full colour), HERMES: 
20% < r{x = 0.023, 5) < 40% (dotted). 



does not contain mixing between different moments. Then, adopting the standard analytical method of full 
moments to the case of the first truncated moment, we have found approximate behaviour of q\ as x — > 0. 
In this way the inverse Mellin transform performed with use of the method of steepest descent implies the 
result q\ IM within our modified "moment of moment" approach. 

We have shown that for small-xo < 0.01 there is a good agreement between the MM results and the reliable 
predictions, obtained with use of the Chebyshev polynomial method. This agreement occurs independently 
either on the shape of the input parametrisation or on the value of Q 2 . It has been also found, that for 
small xo the accuracy of q^ M and q l :MPR - 30 results are similar, being clearly better than in the case of 
gf MPfl - 4 predictions. 

We have presented results concerning the spin structure function g^ s and the truncated at xo = 0.01 
contribution to the Bjorken Sum Rule. It has been found, that the choice of the input parametrisation has a 
large impact on the evaluation of the low-x contribution to the BSR. This contribution can vary from a few 
percents for the flat (~ const) input to tens percents for the steep (~ x~ ' 5 ) one. Recent experimental data 
confirm the rise of the polarised structure functions at small- a;. However, because of the large uncertainties, 
reliable support or contradiction of the theoretical expectations is still out of reach. 
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